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ABSTRACT 

Determination of the rates at which supernovae of Type la (SNe la) occur in 
the early Universe can give signatures of the time spent by the binary progenitor 
systems to reach explosion and of the geometry of the Universe. 

Observations made within the Supernova Cosmology Project are already 
providing the first numbers. Here it is shown that, for any assumed SNe la 
progenitor, SNe la counts up to m R ~ 23 — 26 are useful tests of the SNe la 
progenitor systems and cosmological tracers of a possible non-zero value of 
the cosmological constant, A. The SNe la counts at high redshifts compare 
differently with those at lower redshifts depending on the cosmological model. 
Flat f^-dominated universes would show a more significant increase of the SNe 
la counts at z ~ 1 than a flat, VL M = 1 universe. Here we consider three sorts 
of universes: a flat universe with H = 65 km s _1 Mpc^ 1 , Vt M = 1.0, = 0.0; 
an open universe with H = 65 km s" 1 Mpc" 1 , VLm = 0.3, I^a = 0.0; and a flat, 
A-dominated universe with H = 65 km s' 1 Mpc^ 1 , Vt M = 0.3, fi A = 0.7). On 
the other hand, the SNe la counts from one class of binary progenitors (double 
degenerate systems) should not increase steeply in the z = to z = 1 range, 
contrary to what should be seen for other binary progenitors. A measurement 
of the SNe la counts up to z ~ 1 is within reach of ongoing SNe la searches at 
high redshifts. 

Subject headings: cosmology: general — supernovae: general 
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1. Introduction 

Supernovae are bright stellar explosions which can now be observed up to redshifts 
z ~ 1. Several programmes (Perlmutter et al. 1997, 1998; Schmidt et al. 1997; Garnavich 
et al. 1997) are presently devoted to discover high-redshift supernovae for cosmological 
purposes. The goal is to determine the geometry of the Universe through the effect that 
Hm, have in the magnitude-redshift relationship of Type la supernovae (SNe la). The 
same programmes provide rates of explosion up to large redshifts through proper evaluation 
of control time and efficiency of detection. Indeed, the Supernova Cosmology Project has 
already achieved the first results after evaluating the data collected in its 1995 and 1996 
discovery runs (Pain et al. 1996, 1997). 

In this Letter we outline the theoretical background for the prediction of those number 
counts for different geometries of the Universe and for various SNe la progenitors. 

The pace at which stars formed in the past (the SFR) and the evolutionary clock to 
explosion of those binary supernovae is reflected in the counts. The evolutionary clock 
spans, in the case of SNe la — as opposed to SNe II — , an important fraction of the age of 
the Universe. 

SNe la counts extending up to apparent red magnitudes m R ~ 23 — 25 could provide 
a good test of the SNe la binary progenitor systems. We also find that SNe la counts 
should be sensitive to a A-dominancy in our Universe through a larger increase at z ~ 1. 
Refinements in the determination of the global SFR and continued SNe la searches at large 
redshifts should furnish, in the near future, an accurate enough basis for both tests. 
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2. Modeling 

The global SFR recently derived for high redshift galaxies in the Hubble Deep Field 
(HDF) by Madau et al. (1996) has been extended to lower redshifts (Madau 1997) in 
accordance with the results of Lilly et al. (1995, 1996). The variation of the SFR with z is 
similar to that of the space density of quasars (Shaver et al. 1996), peaking at z ~ 2, the 
latter showing, however, a steeper slope both prior and after maximum. The above SFR(z) 
at z > 2 has been obtained from the UV luminosity density along redshift. Rowan-Robinson 
et al. (1997) have also derived a SFR(z) from ISO infrared observations of galaxies in the 
HDF. The general shape agrees well with Madau's (1996) data but the inferred values are 
higher, maybe implying that about 2/3 of star formation at high z has taken place shrouded 
by dust. In a different approach, Pei & Fall (1995) have dealt with the determination of 
the global star formation history by tracking the evolution of the global HI contents of the 
Universe as measured from Lya QSO absorption-line systems. The results of Madau et 
al. (1996) agree with those from that very different approach. On the other hand, Lilly et 
al. (1996) provide, from the Canada-France-Hawaii Redshift Survey, an estimate of the 
comoving luminosity density of the Universe over the redshift range < z < 1, which can 
be interpreted in terms of a SFR(z) which agrees with previous estimates of a significant 
increase with z up to z ~ 1.5 — 2. Prospects to further explore the global SFR include 
the number counts of SNe II discovered at high redshift, which should trace the activity of 
star formation in the Universe along z. Such numbers on SN II are obtainable in current 
high-z SNe searches. Given the continuous improvement in the determination of the global 
SFR, we base the present calculation of SNe la number counts on the most recent empirical 
results. We evaluate the rate of explosion of SNe la by convolving their time to explosion 
with the SFR. The rate can then be calculated as: 
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(1) 



where R(t) is the SNe la rate after an instantaneous outburst, and t and r are in the SN 
rest frame. Depending on the progenitor systems, in some cases we would have SNe la with 



only a few 10 7 yr after star formation already, thus becoming the brightest optical events 
at z » 1. SFR(t) is derived from the global SFR(z) (Madau et al. 1996; Madau 1997). 
The adopted SFR(z) (for H = 50 km s^ 1 Mpc" 1 , Vt M = 1.0, Ql^ = 0) is shown in the top 
panel of Figure 1, along with the measured data points and their error bars. Note that 
the data points from Connolly et al. (1997) for the 1 < z <,2 range fall somewhat below 
whereas those from Rowan-Robinson et al. (1997) for the 0.6 <, z < 1 range are significantly 
above the adopted curve. The SFR(z) is transformed according to the geometries of the 
model universes considered. Here we restrict this presentation to three favored models of 
the Universe: Model A, a flat universe with H Q = 65 km s _1 , Qm — 1, — 0; Model B, 
an open universe with H = 65 km s^ 1 Mpc^ 1 , Qm = 0.3, Q\ = 0, and Model C, a flat 
universe with H Q = 65 km s _1 Mpc' 1 , Qm = 0.3, Q\ = 0.7. 

The SFR(z) derived for each case is next transformed into SFR(t) (in the comoving 
frame). The derived global SFR(t), when time-integrated over the whole history of the 
corresponding universe (model A: t = 10.0 Gyr; model B: t = 12.2 Gyr; model C: 
t = 13.6 Gyr), produces the observed stars today (Guzman et al. 1997). 

The global SNe la rates, r SNe i a {t) (yr^ 1 Mpc~ 3 ), for each model universe and family 
of SNe la progenitor systems, are calculated in the comoving frame according to (1), and 
integrated over comoving volume to obtain the expected SNe la counts {yr^ 1 sq.deg~ v ) as a 
function of z. We integrate over dV: 



a peak rate of explosion at 10 9 



yr. In other cases, there can be some SNe la exploding even 
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dV = 



(i + n k H$d 2 M y/* 



d(d M )dQ 



(2) 



where du is the proper motion distance and fife = 1 — VLm ~ is the fractional contribution 
of the curvature to the expansion (Carroll, Press, & Turner 1992). du is related to the 
luminosity distance d^ through d^ = (1 + z) du (Weinberg 1972). rf L is calculated as: 



where sinn stands for sinh if fife > and for sin if fife < (both sinn and the fife terms 
disappear from (3) if fife = 0, leaving only the integral times (1 + z)/H ). 

The dependence on cosmological parameters of the comoving volume derivative dV/dzdVL 
differs from one model of Universe to another. A different cosmological effect is the 
age-redshift (t(z)) relationship for different model Universes, which changes the z at which 
the SNe la rates peak. 

The results are tightly linked to the reliability of the global SFR. An estimate of its 
current uncertainty can be obtained by looking at the data points in the top panel of Figure 
1. Increasingly accurate SFRs should steadily improve the SNe la rate predictions. There 
is, as well, increasing evidence of universality in the slope of the initial mass function (IMF) 
from tests in different metallicity and age environments. Such a "universal" mass function 
would be well represented by the Salpeter (1955) power law with x = 0.86 ± 0.23. Both the 
overall trend to convergence in the estimate of the global SFR(z) and the almost constancy 
of the IMF favor the possibility of deriving global values for the SNe la explosion rates. 

Stellar clock: binary progenitors. All evolutionary models for SNe la progenitors involve 
the accretion of mass by a C+O WD from a companion in a close binary system. Current 
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research on SNe la has discarded some formerly proposed SNe la scenarios. The two 
classes of binary systems considered here as the most likely systems giving SNe la are: 
(1) the merging of two white dwarfs — also called double degenerate scenario (DD) — , 
and (2) the accretion by the WD of H from a less evolved star in a close binary system 
—a single degenerate scenario (SD) — . We call this second case a cataclysmic-like system 
(CLS). The time evolution of the SNe la rates after star formation differs from one family 
of progenitor systems to another, and there are variations within each family. We have 
modeled the SNe la rates by means of the same Monte Carlo code as in previous works 
(Ruiz-Lapuente, Burkert, & Canal 1995,1997; Canal, Ruiz-Lapuente, & Burkert 1996, 
hereafter RCB95,97 and CRB96), and adopting different prescriptions for each evolutionary 
path and for the initial binary parameters. We have also compared our predictions with 
those of other authors: (1) The physical input for the merging of two C+O white dwarfs 
(WDs) is modeled both as in Iben & Tutukov (1984) and in RBC95. We allow for a number 
of different physical descriptions of the binary evolution, and we try different values of the 
common envelope parameter a (a measure of the efficiency with which orbital energy is 
used in envelope ejection). (2) The physical modeling of the explosion of a C+O WD when 
it reaches the Chandrasekhar mass by accretion and burning of H from a main-sequence, 
subgiant, or giant companion is modeled as in the early work by Iben & Tutukov (1984), as 
in CRB96, and finally including the most recently proposed variation of the binary evolution 
by Hachisu, Kato, & Nomoto (1996), who find a "wind" solution for fast hydrogen accretion 
by a WD from a red-(sub)giant companion, but we apply the same evolutionary constraints 
as Yungelson et al. (1996). Particularly relevant to the modeling are the slope of the IMF, 
the distribution of mass ratios q of the secondary to the primary, and the distribution of 
initial separations A in the progenitor binary sytems, together with the prescriptions for 
mass transfer. By trying various physical approaches adopted by different modelers and 
exploring different values of the initial binary parameters (IMF, q, A distributions) and 
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assumptions as to the mass transfer rates, we obtain the characteristic behaviour of each 
stellar clock and the uncertainties in the absolute scale of the rates. 

In Figure 1, a comparison of R(t) in (1) with those from other authors (middle panels: 
different IMF, q and A distributions; bottom panel: different assumptions as to the allowed 
mass transfer rates) is made. From such comparison we concluded in our previous work 
that the evolution of the SNe la rates (rise, peak, decline) along cosmic time, for a given 
class of systems, have broad common features shared by the predictions from different 
authors (Iben & Tutukov 1984; Tutukov & Yungelson 1994; Yungelson et al. 1996; our own 
modeling), which basically describe the clock of the corresponding systems. The Hachisu, 
Kato & Nomoto (1996) solution, here called CLS(W), would allow to have some relatively 
young SNe la in the CLS scenario, and the SNe la rate would increase fast with redshift. 
On the other hand, the DD scenario predicts a flatter increase of the rate from z ~ to 
z~ 1, in all model universes. In all cases, the R(t) shown here would give, for our own 
Galaxy and the present time, SNe la rates which are within the current error bars of the 
observational estimates. 

Versions CLS and CLS(W) in the bottom panel of Figure 1 illustrate the range 
of uncertainty in the time of start of the explosions in the SD scenarios. For the DD 
progenitors we have chosen the most favorable physical assumptions, those that enhance 
the numbers of SNe la at high z, in order to best reproduce the first observational results. 
Any other choices would give lower numbers. In fact, the observed numbers seem to be 
above the predictions for this type of binary progenitor system (but one should stress 
here the uncertainty in the global SFR). A clear feature of the DD scenario is not only 
that the SNe la rates do not increase fast towards higher redshifts but also that they 
are lower than those predicted for the SD progenitors. We should note here that we are 
assuming the universality of the distributions of binary parameters determined for the solar 
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neighbourhood (Duquennoy & Mayor 1991). We have explored, however, the effects that 
plausible changes in the initial q and A distributions would have in the outcome and found 
them to be only moderate (see below). 

Counts as a function of magnitude. We compare the model counts with the observations 
by deriving the dependence of the number counts on m R , the apparent magnitude of the 
SNe la in the R photometric band. In fact, in this transformation the intrinsic dispersion 
in brightness of the SNe la is not a key factor, since we are using the N — mag relationship 
measured in intervals of 0.2 to 0.5 mag for our tests (we could also use N — z, thus avoiding 
the transformation to magnitude). The intrinsic dispersion in magnitude is, in contrast, 
very important for cosmological tests using m(z). 

We assume the SNe la to have an average absolute blue magnitude: 

M B = -18.52 + 5 log(h/0.S5) (4) 

where h = H /100 (Perlmutter et al. 1997). 

The distance modulus for each z is calculated from the luminosity distance d^. The 
apparent m R at maximum as a function of z is determined from: 

m R = M B + 5 log(d L (z)) - 5 + K corr (5) 

where K corr is taken from Kim, Goobar, & Perlmutter (1996) (it includes the full 
transformation from B into R magnitudes). We finally calculate the variation of the SNe la 
rate (yr" 1 sq.deg~ l Amag~ x ) with apparent red magnitude m R (we take Am R = 0.5 mag). 

The results are displayed in the four panels of Figure 2. Shown in the Figure are the 
data points currently available (Pain et al. 1996; Pain, private communication). The dips 
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in the curves are entirely due to the shape of the i^-corrections. The different slopes of the 
curves from one family of progenitors to another mainly reflect how fast the SNe la rate 
declines after reaching its peak value plus the time elapsed between star formation and the 
peak (steeper decline and longer time delay for the CLS and CLS(W) models than for the 
DD model: see the middle and bottom panels of Figure 1), and the absolute values of the 
rates at low redshift (m# ~ 20 — 21) are sensitive to the ages of the corresponding model 
universes (to(A) < t (B) < t (C)). The slopes of the rates when comparing different models 
of the Universe are sensitive to the contribution of a non-zero A. The faster increase of the 
rates with magnitude (redshift) along the model sequence A-C corresponds to the increasing 
comoving volume derivatives (dV/dzdQ) in the respective model universes. Such derivative 
is large if Ha gives a significant contribution to f2. As a check of the sensitivity of the results 
to the choices of the binary model parameters, we have also calculated the dN/dmR-mji 
relationship, for the DD progenitor, adopting different q and Aq distributions (dotted line 
in the middle panel of Figure 1): the final curve is not very sensitive to the various choices 
of the distributions explored in this work (the model is shown by the continuous line in 
the four panels of Figure 2). Taking at their face values the SFR and parameters for the 
SNe la progenitor evolution adopted, and also the two data points, the results for both 
the CLS and CLS(W) systems give better fits to the data than those for the DD systems, 
model universe A being most favored and model C giving the worst fit for any kind of 
system. But those conclusions are preliminary and we must await the reduction of the 
uncertainties in the observed counts and in the global SFR for the intermediate-mass stars 
leading to SNe la. As we see from Figure 2, at present the uncertainty in the SNe la counts 
at m R ~ 22 is of the order of a factor of 2. That is the same as the difference between the 
prediction for the DD model and the lowest one for the SD model, and also similar to the 
range covered by the two extreme predictions (CLS and CLS(W)) for the SD model. Any 
increase in the SFR would almost homologously shift upwards all the count predictions. 



- 11 - 



At higher magnitudes, the range of predictions for the SD models becomes narrower while 
the differences with the DD model increase. Thus, if we knew both the SNe la counts at 
m R ~ 23 and the global SFR to better than a factor 1.5, we could already discriminate 
between the SD and the DD models. The sensitivity to the cosmological parameters VL M 
and fl\ is still low at those magnitudes. If we had determined, for instance, the DD model 
to be the right one, then the SNe la counts at m R = 24.5 (corresponding to z ~ 1) in 
the C universe (Qm = 0.3, ^a = 0.7) should be twice those in the A universe (Qm = 1, 
= 0), the factor of increase from m R = 22 to m R = 24.5 being also twice as large in 
case C as compared with case A (both contrasts would be sharper if the SD model were the 
right one). The high-z SN searches to measure Qm, from the variation of the apparent 
magnitudes of SNe la are now almost reaching z — 1 (Perlmutter et al. 1998; Garnavich et 
al. 1997). From samples large enough to derive the SNe la counts at z ~ 1, the constraints 
on both parameters should be tighter than those obtainable from the counts alone, but the 
latter could provide a supplementary test, since it is based on a different approach: the use 
of a volume effect instead of a redshift-magnitude effect. Counts extending to even higher 
magnitudes would more clearly reveal the geometry of the Universe. Suggested SN searches 
up to infrared magnitudes K ~ 26 — 27 (Miralda-Escude & Rees 1997) would extend well 
beyond that point. 

3. Conclusions and future prospects 

The theoretical results presented here are intended to show the potential of the 
comparison between model predictions and the results of undergoing and future SNe 
la searches. The slopes of the dN/dm R -m R curves are mainly sensitive to the general 
characteristics of each SNe la model while the absolute values depend more on model 
parameters and on the SFR. Hence the interest of data covering a broader m R (or z) 
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range than those currently available. On the other hand, the differences in the predictions 
from any SNe la progenitor assumed, for f^-dominated universes as compared with 
matter-dominated and open universes, become significant at large enough values of tur. 
In a Universe dominated by ^a the number counts of SNe la show a larger increase at 
z = 1 due to the large volume encompassed at that redshift. Useful information would be 
obtained by the measurement of the evolution of the SNe la counts up to z = 1 and beyond 
(possibly in the K-band for higher redshifts). 

The first data obtained by Pain et al. (1996) gave a value of 34.4+162 SNe la 
yr^ 1 sq.deg.~ l in a magnitude range of 21.3 < R < 22.3. These first results were obtained 
from three SNe discovered at redshifts 0.374, 0.420, and 0.354 (SN 1994H, SN 1994al & SN 
1994F) and the small number statistics dominates this very first measurement. A larger 
bulk of data, already obtained, will now reduce the statistical uncertainties, and prospects 
to extend the observations up to higher z are on the way. It is thus too soon to extract 
firm conclusions from the data. Our purpose here is to propose a new useful test to be 
completed in the near future. 
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Fig. 1. — Top panel: the adopted global SFR(z), from Madau (1997), together with different 
observational determinations: triangles from Lilly et al. (1995, 1996), diamonds from 
Connolly et al. (1997), squares from Madau et al. (1996), crosses from Rowan-Robinson et 
al. (1997). Middle panel: the SNe la rates after an instantaneous outburst of star formation, 
for the double-degenerate (DD) model. The continuous line is our own result and the dotted 
line that of Tutukov & Yungelson (1994), who adopt a different IMF and also different q 
and A distributions. Bottom panel: same, for the cataclysmic-like system (CLS) model. 
Continuous line is our own result and dotted line the "wind" solution of Hachisu, Kato, & 
Nomoto (1996), with evolutionary prescriptions from Yungelson et al. (1996) and ourselves 
(see text). 
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Fig. 2. — Variation of the SNe la rates with limiting apparent red magnitude m#, for model 
universes A (top left panel), binned to intervals of 1 mag to compare with the results of 
Pain et al (1996, 1997); model universe A but without the binning (top right panel); model 
universe B (bottom left panel) and C (bottom right panel), for the different SNe la progenitor 
systems considered. Continuous lines correspond to the double-degenerate (DD) systems, 
dotted lines to the cataclysmic-like systems (CLS), and long-dashed lines to the CLS systems 
also, for the "wind" solution of Hachisu, Kato, & Nomoto (1996) (see text). 



